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We present the full 0{a) corrections to single Higgs production in e''"e“ collision. The computation is performed 
with the help of GRACE-LOOP where a generalized non-linear gauge fixing condition is implemented. The numerical 
results are checked by testing their UV and IR finiteness as well as their independence on all five non-linear gauge 
parameters. We find that for a 500 GeV collider and a light Higgs boson of mass 150GeV, the total correction is 
small when the results are expressed in terms of a rather than G^j. For a higher Higgs boson mass of 350GeV, 
the correction is of order —10%. 


1. Introduction 

Single Higgs production will be one of the 
most important processes at the future electron- 
positron linear collider. There are two main 
mechanisms to this reaction: Higgs-strahlung and 
W fusion. At tree-level, the calculation is well 
under control and takes into account the inter¬ 
ference effects between Higgs-strahlung and W 
fusion[l]. The one-loop corrections to the two- 
body ^ ZH Higgs-strahlung process have 
been studied both in the standard model[2] and in 
the minimal supersymmetric standard model [3]. 
However, despite the fact that for linear collider 
energies W fusion becomes the dominant Higgs 
production mechanism a full 0(a) calculation has 
been missing. This is indispensable if one wants 
to carry a precision measurement program. The 
loop correction to the HWW vertex has been 
considered [4] and one could argue that it might 
constitute a good approximation for the fusion 
process, but it rests to see how well this approx¬ 
imation fares in comparison with the full calcu¬ 
lation. Very recently one-loop radiative correc¬ 
tions to the fusion process have been investigated 
within supersymmetry[5,6] but again by only tak¬ 
ing into account the contribution of the fermions 
and sfermions to the H/hWW vertex. In this 
paper, we report the full 0{a) corrections con¬ 
sisting of virtual and soft corrections as well as 
hard photon radiation, to single Higgs produc¬ 


tion in e'''e , including both the fusion and Higgs- 
strahlung processes in the standard model.^ 

2. GRACE-LOOP system 

The Feynman diagrams and the correspond¬ 
ing matrix elements for e“'"e“ —> lyDH were gen¬ 
erated by GRACE-LOOP[7], an automatic calcula¬ 
tion system. In addition, the radiative process 
e“'"e“ —> vvH^ with a hard photon has been cal¬ 
culated by GRACE [8]. The latter process is also 
necessary to complete the total correction. We 
deal with two sets of diagrams, the full set and the 
production set. In the standard model and within 
the non-linear gauge fixing conditions which will 
be discussed later, one counts 12 tree-level and 
1350 1-loop diagrams (with 98 pentagon graphs). 
This defines the full set. When the diagrams 
including scalar/pseudoscalar-electron couplings 
proportional to the electron mass are removed 
from the full set, we obtain the production set 
which consists of 2 tree diagrams and 249 1-loop 
diagrams (with 15 pentagon graphs). The full 
set has been used to make careful checks of the 
numerical calculation in quadruple precision. Af¬ 
ter checking the system with the full set at some 
random points in phase space, the phase space 
integration of the production set is done in dou- 

^In the workshop, O. Tarasov also presented the same 
topics. 
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ble precision to obtain cross sections and distri¬ 
butions. The production set is sometimes used 
in quadruple precision to confirm the numerical 
stability of the result. 

The GRACE-LOOP system, which has a model 
file that describes all the interaction vertices de¬ 
rived from a particular Lagrangian, can gener¬ 
ate all the necessary Feynman graphs together 
with their codes, so that matrix elements can be 
generated for the calculation of the cross section 
and event generation. For loop processes, the au¬ 
tomatic calculation goes through an intermedi¬ 
ate step performing some symbolic manipulation. 
This involves REDUCE[9] or F0RM[10] for handling 
all the Dirac and tensor algebra in n-dimension 
for all the interference terms between tree-level 
and 1-loop diagrams. The one-loop diagram con¬ 
tribution from each loop graph 5 at a phase space 
point, is defined as 

dag (X 23? • T‘'-“ (1) 

where T*'’®® is the tree-level amplitude summed 
over all tree-diagrams. Tg°°P is the one-loop am¬ 
plitude contribution of a one-loop diagram g. The 
Feynman trick for the propagator is also automat¬ 
ically applied at this stage. This is then passed to 
a module of libraries for the loop integration con¬ 
taining the FF[11] packages as well as the in-house 
reduction formulas. The system together with 
the one-loop renormalization set-up has been de¬ 
scribed in [ 12 ]. 

3. Checks of the system 

The set of input parameters for the calcula¬ 
tion is the following. Throughout this paper, the 
results are expressed in terms of the fine struc¬ 
ture constant in the Thomson limit uqed = 
137.0359895 and the Z mass Mz = 91.1876GeV. 
The on-shell renormalization program uses Myy 
as input parameter, nonetheless the numerical 
value of M\y is derived through Ar[13] Mw 
thus changes as a function of Mh- We take the 
set M„=Md=58MeV, Ms=92MeV, Mc=1.5GeV, 
Mb=4.7GeV, which gives a “perturbative” value 

^We include NLO QCD corrections and two-loop Higgs 
effects. We take cy.s(M'^)=0.11d, together with Gg 
= 1.16639xl0-5GeV'-2. 


of a{Mz) compatible with the current experimen¬ 
tally driven value. We also take Mtop=174GeV. 
With theses values, we derive =80.3767GeV 
for M//=150GeV, and Mw=80.3158GeV for 
M//=350GeV. For the Higgs-strahlung subpro¬ 
cess we require a Z-width. We have taken a con¬ 
stant Z-width, r2=2.4952GeV. Unless otherwise 
stated our results refer to the full e“e''" ^ i^DH, 
summing over all three types of neutrinos with, 
for electron neutrinos, the effect of interference 
between fusion and Higgs-strahlung. 

The results of the calculation are checked by 
performing four kinds of tests at some point in 
phase space. For these tests to be passed one 
works with the full set in quadruple precision. 
One hrst checks the ultraviolet finiteness of the 
results. The regulator constant Cuv = l/£ — 
"fE + logdTT, n = 4 — 2e is kept in the matrix el¬ 
ements. When one varies this parameter CjjVj 
the ultraviolet finiteness test gives a result that 
is stable over 30 digits. Infrared (IR) finiteness is 
checked by introducing a hctitious photon mass 
A, treated in the code as an input parameter. 
The sum of loop and bremsstrahlung contribu¬ 
tions is stable over 23 digits when varying A. The 
third check relates to the independence on the pa¬ 
rameter kc which is a soft photon cut parameter 
that separates soft photon radiation (analytical 
formula) and the hard photon performed by the 
Monte-Carlo integration. Gauge parameter in¬ 
dependence of the result is performed as a last 
check through a set of five gauge hxing parame¬ 
ters. For the latter a generalized non-linear gauge 
fixing condition [14] has been chosen. 


>Cgf = 
1 


(9g — ieceA^ — ig cos 6widZ^)W^^ 


+ 5H + ikxs)x^ 




1 


- — ( 2 ) 


The X represent the Goldstone. We take = 
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Table 1 

Numerical size of S/M^*) for each non-linear gauge parameter. of graphs’ means the number of graphs 
depending on C. Except for S, no 3rd order term appears in this process. 


c 

# of graphs 




a 

149 

- 

10-28 

10-30 

/3 

314 

- 

10-31 

10-23 

6 

1059 

10-20 

10-20 

10-26 

k 

122 

- 

10-23 

10-23 

€ 

132 

- 

10-21 

10-30 


= 1 so that no “longitudinal” term ap¬ 
pears in the gauge propagators. Not only this 
makes the expressions much simpler and avoids 
unnecessary large cancellations, but it also avoids 
the need for high tensor structures in the loop 
integrals. The use of hve parameters is not re¬ 
dundant as often these parameters check comple¬ 
mentary sets of diagrams. We introduce a generic 
notation ( to represent either of d, /3, S, k, e. For 
each C the first check is made while freezing all 
other four parameters to be 0. We have also made 
checks with two parameters non-zero. In principle 
checking for 2 or 3 values of the gauge parameter 
should be convincing enough. We in fact go one 
step further and perform a complete test of gauge 
parameter independence. To achieve this we gen¬ 
erate for each non-linear gauge parameter C, the 
values of the loop correction to the total differ¬ 
ential cross section as well as the contribution of 
each one-loop diagram contribution for the hve 
values C = 0,±1,±2. A rapid look at the struc¬ 
ture of the Feynman rules of the non-linear gauge 
leads one to conclude that for ^ lypH each 
contribution is a polynomial of (at most) third de¬ 
gree in the gauge parameter and thus, that each 
contribution, dag in Eq.(l) may be written as 

dCTs(0 = dcr(°) -h CdcrW (2^^(2) ^ 3 j^( 3 ) (3^ 

For each contribution dug, it is a straightfor¬ 
ward matter, given the values of dug for the hve 
input C = 0, ±1,^2, to reconstruct (Tg°’^’^’^^ For 
each set of parameters we automatically pick up 
all those diagrams that involve an explicit depen¬ 
dence on the gauge parameter. The number of 


diagrams in this set depends on the parameter 
chosen. In some cases a huge number of diagrams 
is involved. For the process at hand this occurs 
with the parameter <5 where over 1000 diagrams 
contribute to the sums, as shown in Table 1, and 
thus almost the entire set of diagrams, is involved 
in the check. 

We then verify that the differential cross sec¬ 
tion is independent of C 

dcr = ^dCTg(C) = ^dCT(°). (4) 

a a 

In order to check this, we introduce the following 
notation for each parameter. 




Maxy(|d4‘’l) 


i = 1,2,3. 


( 5 ) 


As seen from Table 1 agreement within 20 to 
30 digits is observed. (The agreement gets better 
if one artihcially gives the electron mass a higher 
value.) The gauge parameter dependence check 
not only tests the various components of the input 
file (correct Feynman diagrams for example) but 
also the symbolic manipulation part and most im¬ 
portant of all the correctness of all the reduction 
formulae and the proper implementation of all the 
n-point functions. This is quite useful when one 
deals with 5-point functions as the case at hand. 


4. Numerical results 

The results we show here is based on the pro¬ 
duction set with all 3 neutrinos species and in¬ 
clude the hard bremsstrahlung part. Integration 
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Figure 1. cos 9 distributions of Higgs boson at 
-y/s = 500GeV and Muiggs = 350 GeV. The solid 
and dashed lines show 0{a) corrected and tree 
level distributions, respectively. 


over all photon energies and angles is thus per¬ 
formed. The effect of radiative corrections are 
presented in the a-scheme. 

The overall correction to the total cross sec¬ 
tion at the linear collider energy (Y^=500GeV), 
is small for a light Higgs mass of order 150 GeV as 
shown in Table 2. The radiative correction factor 
A is defined as 


^tree 

The magnitude of the correction increases 
steadily for higher Higgs masses. For example, 
for a Higgs mass of 350 GeV at a center-of-mass 
energy of 500 GeV the correction reaches —10%. 
This —10% correction is distributed almost evenly 
over the range of scattering angles of the Higgs 
boson, see Fig.l. 

The distribution of the Higgs boson energy is, 
however, affected by the 0{a) corrections only in 
the region of the resonance of the Z—boson. The 
resonance structure of the tree level is distorted 
in the form of a radiative tail as shown in Fig. 2. 


Figure 2. Energy distributions of Higgs boson at 
y/s = 500GeV and Mniggs = 350 GeV. The solid 
and dashed lines show 0{a) corrected and tree 
level distributions, respectively. 


5. Summary 

GRACE/LOOP has been applied to calculate the 
full 0{a) corrections to single Higgs production in 
the standard model. One-loop amplitudes of the 
full 1350 diagrams were generated by this sys¬ 
tem. The non-linear gauge fixing condition has 
been introduced and shown to be a very power¬ 
ful tool to check the consistency of the full set 
of amplitudes and the system itself. At •\/s=500 
GeV, the correction is around —10% for a 350 
GeV Higgs mass. The radiative tail was observed 
for the energy distribution of the Higgs boson. It 
means that for more realistic predictions, one has 
somehow to include multiple photon emissions. 
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Table 2 


Higgs mass dependence for the tree total cross sections, the full 0{a) corrected one and the radiative 
correction factor at -y/i = 500 GeV. 


Mh 

(GeV) 

^tree 

(fb) 

^O(a) 

(fb) 

A 

(%) 

150 

61.12 

60.99 ±0.07 

-0.2 

200 

37.33 

37.16 ±0.04 

-0.4 

250 

21.17 

20.63 ±0.02 

-2.5 

300 

10.76 

10.30 ±0.01 

-4.2 

350 

4.603 

4.184 ± 0.004 

-9.1 
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